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Abstract 

An efficient method is presented as a means of an approximate, analytic time- 
dependent solution of the Fokker-Planck equation (FPE) for the Langevin model 
subjected to additive and multiplicative noise. We have assumed that the dynamical 
probability distribution function has the same structure as the exact stationary one 
and that its parameters are expressed in terms of first and second moments, whose 
equations of motion are determined by the FPE. Model calculations have shown that 
dynamical distributions in response to applied signal and force calculated by our 
moment method are in good agreement with those obtained by the partial difference 
equation method. As an application of our method, we present the time-dependent 
Fisher information for the inverse-gamma distribution which is realized in the FPE 
including multiplicative noise only. 
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1 INTRODUCTION 



The Langvin model is a very important model to describe the diffusion behavior in non- 
equilibrium systems, and it has been widely applied to various phenomena in physics, 
chemistry and biology. The Langevin model is usually transformed to the Fokker-Planck 
equation (FPE) which deals with the probability distribution function (PDF) of a state 
variable pQ. It is generally not possible to obtain analytic solutions of the second-order 
partial equations. Indeed, exact analytical solutions of the FPE are known for only a few 
cases. In most cases, approximate solutions are obtained by using analytic or numerical 
methods. Typical analytic methods are an appropriate change of variables, eigenfunction 
expansion, perturbation expansion, path integral, Green's function, moment method, and 
the continued-fraction method pQ. When no analytic solutions are available, numerical 
methods such as finite-difference and finite-element methods have been employed. 

For some Langevin models subjected to additive noise only, exact solutions have been 
obtained. For the linear Langevin model, the exact dynamical solution is expressed by the 
Gaussian distribution with time- dependent mean and variance of a state variable. For the 
FPE including a nonlinear diffusion term, some authors have obtained exact dynamical 
solutions [2], [3j 0] . The generalized FPEs in which time dependences are introduced in 
drift and diffusion terms have been investigated [I] -[7]. 

When multiplicative noise is incorporated to the Langevin model, the problem becomes 
much difficult [8]. For the linear Langevin model subjected to additive and multiplicative 
noise, the exact stationary solution is available, and it has been considerably discussed in 
connection with the non-Gaussian PDF in the nonextensive statistics [9]-[T3]. An exact 
dynamical solution for the linear Langevin model subjected to multiplicative noise only is 
obtained in Ref. [T3] although it does not represent the stationary solution. Approximate 
dynamical solutions of the linear and nonlinear FPEs subjected to multiplicative noise 
have been discussed with some sophisticated methods such as the polynomial expansion 
of the logarithmic PDF [16J, the linearizing transformation [17] and the direct quadrature 
method for moment solution [18] . 

Numerical methods are powerful approaches when exact dynamical solutions are not 
available. Analytical solutions are, however, indispensable in some subjects. A typical 
example is a calculation of the time-dependent Fisher information which is expressed by 
the derivatives of the dynamical PDF with respect to its parameters. In a recent paper 
[19] , we calculated the Fisher information in a typical nonextensive system described by 
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the linear Langevin model subjected to additive and multiplicative noise. We developed an 
analytic dynamical approach to the FPE combined with the g-moment method in which 
moments are evaluated over the escort probability distribution [10] . The dynamical PDFs 
calculated by our moment method are shown to be in good agreement with those obtained 
by the partial difference equation method (PDEM) [T^]. By using the calculated time- 
dependent PDF, we discussed the dynamical properties of the Fisher information [19]. 
It is the purpose of the present study to extend such an analytical approach so as to 
be applied to a wide class of Langevin model with the use of the conventional (normal) 
moment method instead of the g-moment method. 

The paper is organized as follows. In Sec. 2, we discuss the adopted Langevin model 
and moment method to obtain the dynamical PDF. The developed method has been 
applied to the three Langevin models. We present some numerical calculations of the 
time-dependent PDF in response to an applied signal and force. Section 3 is devoted to 
conclusion and discussion on the dynamics of Fisher information of the inverse-gamma 
distribution. 



2 METHOD AND RESULT 



2.1 Fokker- Planck equation 

We have adopted the Langevin model subjected to cross-correlated additive (£) and mul- 
tiplicative noise (77) given by 



f = F(x) + G(x) V (t)+Z(t) + I(t). 



Here F(x) and G(x) are arbitrary functions of x, I(t) stands for an external input, and 
r](t) and express zero-mean Gaussian white noises with correlations given by 



(v(t) v(t')) 

mm) 
mm) 



a 2 5(t-f), 
(3 2 6(t-t f ), 
ea(3 5{t - 1') 
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where a and (3 denote the strengths of multiplicative and additive noise, respectively, and 
e the degree of the cross-correlation between the two noise. 
The FPE is expressed by [2)1 EEJ E2] 



Of 



p(x,t) 



d_ 

dx 



F{x) + / 



[a 2 G(x)G\x) + ta(3 G'(x)] 



p(x,t) 



+ Q-^{[a 2 G(x) 2 + 2eaPG(x)+{3 2 }p(x,t)} J (5) 

where G'(x) = dG(x)/dx, and = and 1 in the Ito and Stratonovich representations, 
respectively. Although we have adopted the Langevin model for a single variable in this 
study, it is straightforward to extend it to the coupled Langevin model with the use of 
the mean- field approximation [T9j . 

For I(t) = I, the stationary PDF of p(x) is expressed by [32 

•1, 



\np(x) 



X(x) + Y(x) - 1 1 - | j In (J[a 2 G(x) 2 + 2ea(3G(x) + (3 2 } 



with 



X(x) 
Y(x) 



2 / dx 
2 / dx 



F(x) 



a 2 G(x) 2 + 2eaf3G(x) + (3 2 



(6) 
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_a 2 G(x) 2 + 2ea(3G(x) + (3 2 

2.2 Equations of motion for the moments 

An equation of motion for the nth moment is given by 



d(x n ) 



9pM x«dx, 



dt J dt 

= n f (x ft_1 F(z)) + (x n_1 /(t)) + ^ [a 2 (x n ~ 1 G(x)G'(x)) + ea(5 (x n - l G\x)) 
n(n — 1 



(9) 



+ 



a 2 (x n - 2 G{x) 2 ) + 2eaf3 (x n - 2 G{x)) + (3 2 (x n - 2 ) 



where suitable boundary conditions are adopted. For n — 1, 2, we obtain 



at 

d(x 2 ) 

at 



= (F(x)) + (I(t)) + ^[a 2 (G(x)G'(x)) + ea/3(G'(x))}, 

= 2(xF(x)} + 2(xl(t)} + <f)[a 2 (xG(x)G'(x)} + ea(3(xG'(x)}} 
+ a 2 (G(x) 2 ) + 2ea(3(G{x)) + f3 2 . 



(10) 
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Expanding x as x = fi + 5x and retaining up to 0(((5x) 2 )), we obtain equations of motion 
for the average \x [= (x)} and variance cr 2 [= (x 2 ) — (x) 2 } given by [13] 

^ = /o + M 2 + ^ (a 2 [^i + 3(010, + M? 3 y 2 ] + ea/3(<7i + 3^ 2 )) + /(*), (13) 



do 2 



= 2/ lf r 2 + (0 + 1)(# 2 + 2 M?2 )aV + 2ea/3(0 + l)^ 2 
+ a 2 # 2 + 2ea(3g + p 2 , 



(14) 



where /, = (l/£\)d^F( t i)/dx t and g e = (l/i\)d e G(fi)/dx e . 



2.3 Model A 

2.3.1 Stationary distribution 

Our dynamical moment approach will be applied to the three Langevin models A, B and 
C, which will be separately discussed in Sees. 2,3, 2.4 and 2.5, respectively. 
First we consider the model A in which F(x) and G(x) are given by 

F(x) = -Ax, (15) 
G(x) = x, (16) 

with e = 0.0 (i.e., without the cross-correlation), where A expresses the relaxation rate. 
The model A has been adopted as a microscopic model for nonextensive systems |ll]-|13j. 
From Eq. (J5J), the FPE in the Stratonovich representation is given by 

d d ( (3 2 \ d 2 

dt P ^ = ^ [ Ax ~ J WMM) + (y J Q^Pfat) 



+ 



'a 2 



_ A 
2~j dx 



x-^{xp(x,t)} 



(17) 



v(x) = ,1)53^^1, (18 ) 



By using Eqs. (jSI)-©, we obtain the stationary PDF given by [19 



with 



a 

(2X + a 2 ) 
2c7 2 ' 

/ 

a/3' 

/7Fr(6)r(6-_|) 
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(19) 

(20) 
(21) 

(22) 



a | r(6 + ic 

By using Eq. (fl8l) . we obtain the mean and variance in the stationary state given by 

~ a(b-l) ~ (2\-a 2 Y (23) 
[{b-lf + c 2 } [a 2 !, 2 + (3 2 ) 
a 2 (b- 1) 2 (26 -3) 2(A-a 2 ) ' 1 ' 

Depending on the model parameters, the stationary PDF given by Eq. ([18]) may reproduce 
various PDFs such as the Gaussian, g-Gaussian, Cauchy and inverse-gamma PDFs [19J. 
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2.3.2 Dynamical distribution 

It is worthwhile to remind the dynamical solution of the FPE given by Eq. (ITT)) in the 
limit of a = 0.0 (i.e., additive noise only), for which the time-dependent solution is given 
by 



1 



p ( Xj t ) = - : e -\*-m?iMty ; (25 ) 



2 



ft™®* 

with fj,(t) and cr(t) 2 satisfying equations of motion given by 

dfi(t) 



dt 
da(t) 



2 



Xfi(t) + I(t), (26) 
2Xa(t) 2 + (3 2 . (27) 



dt 

In order to derive the dynamical solution of the FPE for a ^ 0.0 given by Eq. (1171) . 
we adopt the moment approach with the following steps: 

(1) We assume that dynamical PDF has the same structure as the stationary one, as given 
by 

/ J_\ exp[2c(t) tan-HaftM] ^ 
P[X ' t} ~ \Z(t)J [l + a(i) 2 x 2 p) ' 1 j 

with 



^ <t) \T[b(t)+ic{t)} V (29) 

(2) With the assumption (1), we first tried to derive equations of motion for the parameters 
of a(t), b(t) and c(t), by using the FPE after Refs. [21 El II]- Unfortunately, it did not 
work because functional forms in the left and right sides of the FPE become different. 

Then we tried to express the parameters in terms of importance quantities of fj,(t) and 
u(t) 2 such as to be consistent with the relations for the stationary state given by Eqs. 
(123]) and (1241) . Because the number of parameters (three) is larger than two for /i(t) and 
cx(i) 2 , the parameters of a(t), b(t) and c(t) cannot be uniquely expressed in terms of /i(t) 
and cr{t) 2 from Eqs. (!23l and ff24l) . If the first three moments in the stationary state are 
available, it is possible to uniquely express a(t), b(t) and c(t) in terms of them, though 
such a calculation is laborious. 

In order to overcome the above problem, we have imposed an additional condition 
that expressions for the parameters yield the consistent result in the two limiting cases 
of a — > and j3 — > 0. After several tries, we have decided that b(t) and c(t) in Eqs. ff28l) 
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and (129]) are expressed as 

[1 + q 2 {Mt) 2 + 3a(t) 2 }] 
* W " 2aV(t) 2 ' m 

C{t) = 2^W ' ( } 

with the time-independent a (= a/ (3) given by Eq. (1T91) . The relations given by Eqs. (1501) 
and (1311) are consistent with Eqs. (1231) and (124j) for the stationary state and they satisfy 
the above-mentioned limiting conditions, as will be shown shortly 

(3) Equations of motion for fi(t) and <r(t) 2 in Eqs. (13"U1) and (13T1) are obtained from Eqs. 
(TT3"j) - (fTEj) . as given by 

d Jf = - Mt)+m + eMl, (32) 

^2! = -2Aa(() 2 + 2aV(i) 2 + Q'V<) 2 + )3 2 - (33) 

Thus the dynamical solution of the FPE given by Eq. (|T7|) is expressed by Eqs. (128]) - (l33l) . 

In the following, we will show that the relations given by Eqs. (1301) and (13T1) lead to 
results consistent in the two limiting cases of a — > 0.0 and (3 —>■ 0.0. 

(a) a — > case 

In the limit of a — > 0.0 (ie., additive noise only), p(x,t) given by Eq. (1281) reduces to 

pfot) OC e -«W^W- 2 +2a(t)c(t). (34) 

because Eqs. ( 130]) and (l3Tj) with a — > 0.0 yield 



2(x(t) 2 2<r(t) 



2a(t)c(() = [l + aWyOT],(^ ^ (3?) 
Equation (135]) agrees with the Gaussian distribution given by Eq. fl25|) 
(b) p -»• case 

In the opposite limit of /3 = 0.0 (ie., multiplicative noise only), the stationary PDF 
given by Eqs. (TIS1) and (1221) with / > leads to the inverse-gamma distribution expressed 
by 

p(x) = x- s e~^ e(x), (38) 
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where 



5 = 2b, (39) 

k = — . (40) 
a 

Here T(x) denotes the gamma function and 0(t) the Heaviside function: Q(t) = 1 for 
t > and zero otherwise. From Eq. (I3"8l) . we obtain the average and variance in the 
stationary state given by 

" = w^y (41) 

" 2 = (s-2) 2 (s-sy (42) 

from which 5 and k are expressed in terms of /x and a 2 as 

/i 2 + 3a 2 



a 2 



(43) 



/i 2 + cr 2 

« = 2 M- 44 

On the contrary, the dynamical PDF given by Eq. ff28l) in the limit of (3 — > 0.0 (and 
I > 0) reduces to 

g—2c(t) / a(t)x g— K,(t)/x 

p(?> * x2b( t) e ( x ) -> 0(a;) ' (45) 

because Eqs. f )30|) and (131]) with /3 — > 0.0 (a — ► oo) lead to 

xm [i + « 2 {M^) 2 + 3a(t) 2 }] ^{tf + 3tx(t) 2 

5(t) = 2m = — — ' (46) 

K (t) = 2C(^) _ [1 + g!Mg + glffUMg 5 /i(^) 2 + ^ (f y (47) 

a(i) a 2 cr(t) 2 °"(^) 2 

Equations <H5]> . (1461) and (@ZD agree with Eqs. (15511 . (gig) and (jUl) . respectively. 

Thus the expressions given by Eqs. ( J301) and (l31il yield the consistent result covering 
the two limits of a — > 0.0 and (3 —* 0.0. 

2.4 Model B 

2.4.1 Stationary distribution 

Next we consider the model B in which F(x) and G(x) are given by Eqs. ( 1T51) and ( |T6i) 
with e 7^ 0. From Eqs. (jSJ)-©, the stationary PDF is given by 

_ / 1 \ exp[2ctan~ 1 {a(a: + /)}] 
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with 



a 

m^< (49) 

(50) 



2a 2 



af}\Jl — e 



2 



(51) 



/ = 'A (52) 

_ 0Fr(6)r(6-i) 

Z " a | T(6 + ic) i 2 ■ (53) 
Equations (|4"S1) and ( 1531) yield the average and variance in the stationary state given by 

c , (2I + ea(3) 

" = ^T)- / = ^A^T' (54) 

[(ft ~ I) 2 + c 2 ] (<» V + 2ea/?/* + /? 2 ) 
a a 2 (6- l) 2 (26-3) 2(A-a 2 ) ' 1 } 

2.4.2 Dynamical distribution 

With the use of the procedure mentioned for the model A in Sec. 2.3, the dynamical PDF 
p(x, t) of the model B is assumed to be given by Eq. fj48l but with b and c replaced by 

= MH ^f? +M tn (56) 
m ^+a 2 {W)+f] 2 +°m]W)+n 

C[t) ~ 2aa(ty ' 

which agree with Eqs. (|54|) and (|55|) in the stationary state. Equations of motion for n(t) 
and cr(t) 2 in Eqs. ( 1561) and (1571) are given by 

^ = -W)+I(t) + ^ + ^, (58) 
(x^ 2 2 

= -2Aa(t) 2 + 2c?a{tf + aVW 2 + 2ea/fy(t) + (3\ (59) 

which are derived from Eqs. ffTBl - ffTol) . 

It is noted that in the limits of a — > 0.0 and (3 — > 0.0, the cross-correlation between 
additive and multiplicative noise does not work, and the result for the model B reduces to 
that for the model A. Thus the moment method with the use of Eqs. (156|) and (l57j) leads 
to the results consistent in the limits of a — > 0.0 and (3 — > 0.0, where p(x,t) becomes the 
Gaussian and inverse-gamma distributions, respectively. 
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2.5 Model C 

2.5.1 Stationary distribution 

Now we consider the model C in which F(x) and G(x) are given by 

F(x) = -\(x + s), 
G(x) = Vx 2 + 2sx + r 2 , 



(60) 
(61) 



with e = 0.0 where A expresses the relaxation rate, and r and s are parameters. We 
assume that e = 0.0 because we cannot obtain the analytic stationary PDF for e ^ 0.0. 
The model C with r = s = 0.0 is nothing but the model A. From Eqs. (jEJ)-(jEJ), the 
stationary PDF for the model C is given by 



p(x) oc 



1 - —(x + s) 



,Ylx) 



where 



Y(x) 



S cty/D, 
' 21 



In 



x + s - \[D 



x + s + y/D 
1 ( a(x + s 



with 



21 

a 2 (x + s) 



D 



tan 



for D > 
for D < 
for D = 



a 2 (s 2 -r 2 )-(3 2 
(2A + a 2 ) 



2a 2 

When we consider the case of D < 0, the stationary PDF is rewritten as 

1 \ exp[2ctan _1 {a(x + s)}] 



p(x) 



[l + a 2 (x + s) 



with 



a 



f3 2 + a 2 {r 2 - s 2 ) 
I 



a y P 2 + OL 2 {r 2 — s 2 ) 

/?r(6)r(6-|) 



a | r(6 + ic) | 2 
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(62) 

(63) 
(64) 
(65) 

(66) 
(67) 

(68) 

(69) 
(70) 
(71) 



From Eq. (1681) . we obtain \i and a 2 in the stationary state expressed by 

" = l^)- s = W^)'^ (72) 

[(b - l) 2 + c 2 ] [ a y + 2g/ i + r 2 ) + /? 2 ] 
a a 2 (6- l) 2 (26-3) 2(A-a 2 ) ' 1 j 

2.5.2 Dynamical distribution 

In order to obtain the dynamical PDF for the model C, we adopt the same procedure 
as those for the models A and B. We assume that the dynamical solution p(x, t) of the 
model C is given by Eq. fl68|) but with b and c replaced by 

= [1 +n 2 {[Mt) + sl' + 3<7(t) 2 }] ^ 
2a A a\ty 

_ [l + a 2 {[fi(t) + s] 2 + a(t) 2 }Mt) + s] 
C[t) ~ 2^W ' [7b) 

which agree with Eqs. (1721) and (1731) in the stationary state. Equations of motion for fi(t) 
and cr(t) 2 in Eqs. ( 1741 and (175]) are given by 

dj f- = - (a -y)w«) +»]+/(«), (76) 
*f! = _ 2(A _ aV((f + aW + 2s/l(()+r2 ] +/32 , (77) 

which are derived from Eqs. ( fTBl . ( Il4l) . ( |60l) and ( loTT) . It is easy to see that the moment 
method with the use of Eqs. ( 1741) and ( 1751) lead to the results consistent in the limits of a — > 
0.0 and (3 — > 0.0, where p{x,t) becomes the Gaussian and inverse-gamma distributions, 
respectively. 

2.6 Model calculations 

We will present some numerical calculations in this subsection. In order to examine the 
validity of the moment approach, we have employed the partial difference equation derived 
from Eq. ([5]) with <fi — 1, as given by 



p(x, t + v) = 
+ 



P(x, t) + (-F' + y [(C) 2 + GG&] + 6 -f-G^) v p(x, t) 

\2~u) \p(x + u) -p(x-u)} 



—F - I(t) + ^GG> + ^-G' 
w 2 2 



( y G 2 + ea/?G + y ) ) [p(* + «,*) + p(z - u, t) - 2p(x, t)} , (78) 
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where u and v denote incremental steps of x and t, respectively. We impose the boundary 
condition: 



p{x,t) = 0, 



for I x |> x. 



rn 



(79) 



with x m = 5, and the initial condition of p{x, 0) = po(x) where po(x) is the stationary 
PDF. We have chosen parameters of u = 0.05 and v = 0.0001 such as to satisfy the 
condition: (a 2 x 2 n v/2u 2 ) < 1/2, which is required for stable, convergent solutions of the 



with AI = 0.5 and I b = 0.0 to the model B with A = 1.0, a = 0.5, (3 = 0.5 and e = 0.5. 
Figure [1] shows the time-dependence of the PDF at various t in response to an applied 
input. Solid curves express the results of the moment method calculated with the use of 
Eqs. (HHD, ([H3D, d56])-([59]), and dashed curves denote those of the PDEM with Eq. (P751). 
Figure [2] shows the time-dependence of p,(t) and cr(t) 2 calculated by the moment method 
with Eqs. ( 1551) and (1591 (solid curves) and by the PDEM with Eq. (175)) (dashed curves). 
At < t < 2.0 where no input signal is applied, we obtain p,(t) = 0.071 and <j{t) 2 = 0.179. 
The PDF at t < 2.0 is not symmetric with respect to its center because of the introduced 
correlation of e = 0.5. By an applied pulse at 2.0 <t< 6.0, fi(t) and u{t) 2 are increased, 
and the position of p(x, t) moves rightward with slightly distorted shapes. After an input 
pulse diminishes at t > 6.0, the PDF gradually restores to its original stationary shape. 

Next we apply the pulse input given by Eq. (1801) with AI = 0.5 and I b = 0.0 to 
the model C with A = 1.0, a = 0.5, (3 = 0.5, r = 0.5 and s = 0.2. Figure [3] shows 
the time-dependent PDF at various t in response to an applied input. Results of the 
moment method calculated with the use of Eqs. (1651) . (171]) . (173]) - (I77I) are shown by solid 
curves while those of the PDEM with Eq. (I75|) are expressed by dashed curves. The time- 
dependent fi(t) and cr(t) 2 are shown in Fig. HI solid curves denote the results calculated 
by the moment method with Eqs. (1761) and (177j) : dashed curves express those by the 
PDEM with Eq. ( 175|) . In the stationary state at t < 2.0, we obtain fi(t) = —0.2 and 
cr(t) 2 = 0.2. By an applied pulse at 2.0 < t < 6.0, fi(t) and u(t) 2 are increased and the 
position of p(x, t) moves rightward with slightly changed shapes. 

It is possible to calculate the response to temporal changes in the model parameters 
such as A, a and (3. As an example, we introduce the time-dependent relaxation rate \{t) 



PDEM. 



First we apply a pulse input signal given by 



I(t) = AI 6(t - 2)9(6 -t) + I h , 



(80) 
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given by 

X(t) = AA Q(t - 2)6(6 — t) + X b , (81) 

with AA = 1.0 and X b = 1.0 to the model C with a = 0.5, p = 0.0, r = 0.5, s = 0.2 
and / = 0.0. Equation ( 1HTT) stands for an application of an external force of — AA x at 
2.0 < t < 6.0. The time dependent p(x, t) is plotted in Fig. [51 where solid and dashed 
curves denote the results of the moment method and PDEM, respectively. Figure [6] shows 
the time dependences of /i(t) and c{t) 2 . When an external force is applied at 2.0 < t < 6.0, 
the width of the PDF is reduced and cr(t) 2 is decreased while fi(t) has no changes. 

It is noted in Figs. [T]l6]that results of p(x, t), nit) and cx(t) 2 calculated by the moment 
method are in good agreement with those obtained by the PDEM. 

3 CONCLUSION AND DISCUSSION 

The moment approach to the FPE discussed in preceding Sec. 2 may be applied to various 
Langevin models provided analytic expressions for the PDF and for the first- and second- 
order moments in the stationary state are available. For example, when F(x) and G(x) 
are given by 

F(x) = -Xx\x\ r ~\ (82) 
G{x) = x\x\ s ~\ (83) 

for r > 0, s > 0, the stationary PDF with e = 0.0 is given by [T2"| [T3] 

p(x) oc (a 2 |x| 2s + (3 2 )~ 1/2 exp[X(x) + Y(x)}, (84) 

with 

, > /2A|x| r+1 \^/ r + 1 r + 1 a 2 \x\ 2s \ , . 

TW , (2I\x\\ 11 a 2 \x\ 2s \ , , 

where F(a, b, c; z) denotes the hypergeometric function. Equations of motion for fi and 
cr 2 are given by [13] 



/ « 2 ^ 



2 



(yj[3/i|/,| 2s - 2 + S ( S -l)(2 S -l) / ,| / i| 2 -V 2 ], (87) 

-2Ar|/i| r -V 2 + 2s(2s - l)a 2 \n\ 2s - 2 a 2 + a 2 \/j\ 2s + (3 2 . (88) 
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If analytic expressions for stationary values of /i and a 2 are obtainable from Eqs. (jSj 
and ([88]) . we may apply our moment method to the FPE given by Eqs. (jSJ), (1821) and 
(185]) with the following steps: (1) adopting the stationary PDF given by Eqs. (!84]) - (l86l) . 
and (2) expressing its parameters in terms of the time-dependent /i(t) and cr(t) 2 in an 
appropriate way, as mentioned for models A, B and C. 

As an application of our method, we have calculated the Fisher information for the 
dynamical inverse-gamma distribution, which is realized for (3 = 0.0 in the model A [Eq. 

(USD], 

K S(t)-l 

T[5{t) - 1\ 

with the time-dependent S(t) and n(t) given by Eqs. (l4"6l) and (]47j) . With the use of Eq. 
(189]) . the Fisher information matrix given by 



is expressed by 



985 = ^[S(t)-l], (91) 

9KK = K{tf ' (92) 

K{t) 

where ip(x) and ip'(x) are di- and tri-gamma functions, respectively. 

Figure [7] shows the time-dependent inverse-gamma distribution p(x, t) when an input 
pulse given by Eq. (IHOl with AI = 0.3 and h = 0.2 is applied to the model A with A = 1.0, 
a = 0.5 and (3 = 0.0. Solid and dashed curves show the results of the moment method and 
the PDEM, respectively. The time dependences of /i(t) and &{t) 2 and of Fisher information 
are plotted in Fig. [8J(a) and (b), respectively, where solid (dashed) curves express the result 
of the moment method (PDEM). By an applied pulse at 2.0 < t < 6.0, fi(t) and u{t) 2 
are increased and the position of the PDF moves rightward with an increased width. An 
applied pulse increases g KK while it decreases g& K . An interesting behavior is observed 
in gss which is decreased at t = 2.0 when a pulse is applied, but afterward it seems to 
gradually reduce to the stationary value. A similar behavior is realized in g$$ also at 
t > 6.0 when the applied pulse is off. 

In our previous paper [19], we applied the q- moment approach to the model A, deriving 
equations similar to Eqs. (I30p - (|55]) in the (normal) moment approach. There are some 
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differences between the q- and normal-moment approaches. The stationary variance a 2 in 
the g-moment approach evaluated over the escort distribution is stable for < a 2 / A < oo 
whereas a 2 in the normal-moment approach is stable for < a 2 /X < 1.0 [Eqs. (!24j) . (1551) 
or (175]) ]. Although the time dependences of (i g (t) and cr(t) 2 calculated by the g-moment 
approach are similar to those of and u(t) 2 for a small a, the difference between 
them becomes significant for a large a (see Fig. 13 in Ref. [IH]). These differences yield 
the quantitative difference in p(x, t) calculated by the normal- and g-moment methods, 
although both the methods lead to qualitatively similar results. We note that in the limit 
of a = 0.0 (i.e., additive noise only), the dynamical solution given by the q- or normal- 
moment method reduces to the Gaussian solution given by Eqs. (!25l) - (!27j) . Thus the q- 
or normal-moment approach is a generalization of the Gaussian solution to the FPE with 
a 7^ 0.0 given by Eq. (j^D- 

In summary, by using the second-order moment method, we have discussed the analytic 
time-dependent solution of the FPE which includes additive and multiplicative noise as 
well as external perturbations. It has been demonstrated that dynamical PDFs calculated 
by the moment approach are in good agreement with those obtained by the PDEM. Our 
moment method has some disadvantages. The variance a 2 diverges at a 2 /X > 1.0, for 
which our method cannot be applied. If an applied perturbation induces a large fi(t) 
and/or cr(t) 2 , our method leads to poor results which are not in good agreement to those 
calculated by the PDEM. These are inherent in the moment approximation in which 
each moment is required to be small. Despite these disadvantages, however, our moment 
method has following advantages: (i) obtained dynamical solutions are compatible with 
the exact stationary solutions in the Langevin models A, B and C, (ii) it is useful for 
various subjects in which analytical dynamical PDFs are indispensable (e.g., Ref. [19]), 
and (iii) the second-order moment approach is more tractable than sophisticated methods 
[IS]- [IB] for the FPE subjected to multiplicative noise. As for the item (iii), it is possible 
to take account of contributions from higher-order moments than second-order ones with 
the use of Eq. ({TO]) , though actual calculations become tedious. 
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Figure 1: (Color online) The time-dependent probability distribution p(x, t) of the model 
B with A = 1.0, a = 0.5, (3 = 0.5 and e = 0.5, in response to an applied pulse given by 
Eq. ( ISDj) with AI = 0.5 and = 0.0: solid and dashed curves denote results calculated 
by the moment method and the PDEM, respectively, curves being consecutively shifted 
downward by 0.25 for a clarity of the figure. 



Figure 2: (Color online) The time dependence of /i(t) and u(t) 2 of the model B with 
A = 1.0, a = 0.5, (3 = 0.5 and e = 0.5, calculated by the moment method with Eqs. (1761) 
and (177)) (solid curves) and by the PDEM (dashed curves) in response to an applied pulse 
I(t) given by Eq. (1801) with AI = 0.5 and = 0.0 (the chain curve). 



Figure 3: (Color online) The time-dependent probability distribution p(x, t) of the model 
C with A = 1.0, a = 0.5, (3 = 0.5, r = 0.5 and s = 0.2, in response to an applied pulse 
given by Eq. (180)) with AI = 0.5 and lb = 0.0: solid and dashed curves denote results 
calculated by the moment method and the PDEM, respectively, curves being consecutively 
shifted downward by 0.25 for a clarity of the figure. 



Figure 4: (Color online) The time dependence of /i(t) and u{t) 2 of the model C with 
A = 1.0, a = 0.5, (3 = 0.5, r = 0.5 and s = 0.2, calculated by the moment method with 
Eqs. (1751) and (1771) (solid curves) and by the PDEM (dashed curves), in response to an 
applied pulse I(t) given by Eq. (IBOl) with AI = 0.5 and = 0.0 (the chain curve). 



Figure 5: (Color online) The time-dependent probability distribution p(x, t) of the model 
C with I = 0.0, a = 0.5, (3 = 0.0, r = 0.5 and s = 0.2, when the time-dependent relaxation 
rate \(t) given by Eq. (IHT!) is applied: solid and dashed curves denote results calculated 
by the moment method and the PDEM, respectively, curves being consecutively shifted 
downward by 0.25 for a clarity of the figure. 



Figure 6: (Color online) The time dependence of /i(t) and u(t) 2 of the model C with 
A = 1.0, a = 0.5, [3 = 0.5, r = 0.5 and s = 0.2, calculated by the moment method 
with Eqs. (176)) and (ITT}) (solid curves) and by the PDEM (dashed curves) when the time- 
dependent relaxation rate \(t) given by Eq. (IHT)) is applied (the chain curve), A being 
divided by a factor of five. 
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Figure 7: (Color online) The time-dependence of the inverse-gamma distribution p(x, t) 
of the model A with A = 1.0, a = 0.5 and ft = 0.0 in response to an applied pulse given by 
Eq. flHOj) with AI = 0.3 and lb = 0.2: solid and dashed curves denote results calculated 
by the moment method and the PDEM, respectively. Curves are consecutively shifted 
downward by 1.0 for a clarity of the figure. 



Figure 8: (Color online) (a) The time dependence of /i(t) and a{t) 2 of the model A 
with A = 1.0, a = 0.5 and ft = 0.0 calculated by the moment method (solid curves) and 
PDEM (dashed curves) in response to an applied pulse given by Eq. (18*01) with AI = 0.3 
and h = 0.2 (the chain curve), u{t) 2 being multiplied by a factor of five, (b) The time 
dependence of the Fisher information matrix calculated by Eqs. fl9Tl) - fl93l : gss (the solid 
curve), g KK (the dashed curve) and g$ K (the chain curve), gss and g KK being multiplied by 
factors of 10 and 1/10, respectively (see text). 
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